#ifndef DEFINE_H
#define DEFINE_H

// #define A[i][j] A[i*Nx + j]
// #define A_PLUS[i][j] A[(i+1)*Nx + j]
// #define A_MINUS[i][j] A[(i-1)*Nx + j]
// #define _A[i][j] A[(i)*Nx + j + 1]
// #define _A_PLUS[i][j] A[(i+1)*Nx + j + 1]
// #define _A_MINUS[i][j] A[(i-1)*Nx + j + 1]

#define N 1000
#define Nx (N+1)
#define Nz 2000
#define L 1.
#define Lz 1.
#define R 0.1
#define Y 1.
#define D 1.
#define X0I (N/2)
#define X0 ((double) X0I/N)
#define ITMAX 100
#define EPSz 0.000000001
#define EPSx 0.02
#endif